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Abstract 

Let (v,p) be a smooth solution pair of the velocity and the pressure 
for the Navier-Stokes(Euler) equations on R N x (0, T), N > 3. We set 
the Bernoulli function Q = ^\v\ 2 +p. Under suitable decay conditions 
at infinity for (v,p) we prove that for almost all a(t) and f3(t) defined 
on (0, T) there holds 

/ (\7rM 2 + v\u\ 2 \ dx = v I \VQ\dS 

J{a(t)<Q(x,t)<p(t)} \ 2(Jt / J{Q(x,t)=0(t)} 

-v I \VQ\dS, 

J{Q(x,t)=a{t)} 

where u = curl v is the vorticity. This shows that, in each region 
squeezed between two levels of the Bernoulli function, besides the en- 
ergy dissipation due to the enstrophy, the energy flows into the region 
through the level hypersurface having the higher level, and the energy 
flows out of the region through the level hypersurface with the lower 
level. Passing a(t) | inf^gjgjv Q(x, t) and f3(t) t sup^g^jv Q(x, t), we re- 
cover the well-known energy equality, J* RJV M 2 = — v J r n |w| 2 clx. A 
weaker version of the above equality under the weaker decay assump- 
tion of the solution at spatial infinity is also derived. The stationary 
version of the equality implies the previous Liouville type results on 
the Navier-Stokes equations. 



1 



AMS Subject Classification Number: 35Q30, 35Q31,76D05 

Key words: Navier-Stokes equations, Euler equations, localized energy 

equality 



1 Introduction 

We are concerned here on the incompressible Navier-Stokes(Euler) equations 
on R N , N > 2. 



where v = v(x,t) = (vi(x, t), • • • ,VN(x,t)) is the velocity, and p = p(x,t) 
is the pressure. We assume the viscosity satisfies v > 0. In the case v > 
the system (NS,E) becomes the Navier-Stokes equations, while for v = the 
system (NS,E) is the Euler equations. By the system (NS,E) we represent 
both of the cases of the Navier-Stokes and the Euler equations. Below we 
denote the vorticity of the vector field v in M. N defined by u = {djVk — 
dkVj}f k=1 .j >k , the magnitude of which is given by 



As is well-known in most of the fluid mechanics text books, the smooth 
solution v of the system (NS,E) with sufficiently fast decays at spatial infinity 
satisfies the energy equality, 



Our aim in this paper is to localize the domain of the integration in (11. ip . 
The domains are characterized by their boundaries, which are the level hyper- 
surfaces of the Bernoulli function, 



Our localized integral equalities refine the classical equality ( 11. ip . in the sense 
that under suitable integrability conditions for the solutions and particular 






(1.1) 



Q(x,t) := -\v (x,t)| 2 +p{x,t). 



(1.2) 
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choice of the levels of the Bernoulli function we recover ( II. ip . The first the- 
orem below concerns on these localized energy equalities under milder con- 
ditions on the asymptotic behavior for the solution (v,p) at spatial infinity. 



Theorem 1.1 Let N > 2, and (v,p) be a smooth solution of (NS,E) on 
M. N x (0, T). Let Q(x,t) be defined as in M.2\) . Suppose there exists Qo = 
Qo(t) such that 

lim Q(x,t) = Q (t) (1.3) 



uni 



iformly for each t G (0, T). Then, for almost every a(t),/3(t) defined on 
(0, T) such that either j3(t) > a(t) > or > (3{t) > a(t) we have the 
following equality. 



L 



{a(t)<Q(x,t)-Q (t)<l3(t)} 



-— !) + v\oj\ ax 



2dt 

{Q(x,t)-Q (t)=/3(t)} J {Q(x,t)-Q (t)=a(t)} 



\VQ\dS-v [ \VQ\dS. 

(1.4) 



Remark 1.1 Physically the equality (ll.4p shows new energy balance in each 
region squeezed between two levels of the Bernoulli function. Besides the 
well-known energy dissipation due to the enstrophy the energy flows into the 
region through the level hypersurface having the higher level, and the energy 
flows out of the region through the level hypersurface with the lower level. 
This tendency of flow of the energy into lower values of the Bernoulli func- 
tion is proportional to the viscosity and the magnitude of the gradient of the 
Bernoulli function, and there exists no such phenomena in the inviscid case 
of the Euler equations, where energy is just conserved on the average in each 
strip between the two levels of the Bernoulli function. 

Remark 1.2 If we approach f3(t) t suPxgiR^ Q{ x ,t) ~ Qo(t), then, since 
f{Q( x t)-Qo{t)=p(t)} 0, we obtain from ( II. 4p that 



' d v\ 2 + u\uj\ 2 \dx = -v I \VQ\dS < 0. 



{a(t)+Q (t)<Q(x,t)} \2 9t / J {Q(x,t)=a(t)+Q (t)} 



(1.5) 
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Similarly, if we approach a(t) j. inf^g^iv Q(x, t) — Qo(t), then we obtain from 
( OP that 



' d 'v\ 2 + v\uj\' 2 ) dx = v I \VQ\dS>0. 



{Q(x,t)<0(t)+Q o (t)} \2 9t / J {Q(x,t)=P(t)+Q (t)} 

(1-6) 

Remark 1.3 If we let a(t) and (3(t) approach respectively to the criti- 
cal values of Q(x,t) — Qo(t) with a(t)(3(t) > 0, then the righthand side 
of (II. 4p vanishes, and the energy balance holds in the region of the set 
{iGM w | a(t) < Q(x, t) - Q Q (t) < f3(t)}. 

The following theorem derives stronger version of equalities than (jl.4p in 
the sense that there is no sign condition on the functions a (£),/?(£), but un- 
der the stronger assumptions on the behaviors of the solutions at spatial 
infinity represented by the integrability of the solutions. 

Theorem 1.2 Let N > 3, and (v,p) be a smooth solution of (NS,E) on 
M. N x (0, T). Suppose (v,p) satisfies the following conditions: there exists 
Po(t) defined on (0, T) such that 

lim (\v(x,t)\ + \p(x,t)- Po (t)\) = 0, (1.7) 

|x|— >-oo 

uniformly for each t G (0, T), 

v E L°°(0,T;L^(R N )nL^(R N )), (1.8) 

and 

^|f| 2 G L°°(0, T; L 1 (IR Af )). (1.9) 

Then, for almost all real valued functions a(t) < (3(t) defined on (0,T) we 
have the following equalities. 



I, 



{a(t)+po(t)<Q(x,t)<P(t)+p (t)} 



1 9 , l2 ir'l, 

~-\v\ +„M )dx 



v / \VQ\dS-v / \VQ\dS. 

1 {Q(x,t)=0(t)+p o (t)} J {Q{x,t)=a(t)+p (t)} 



1.10) 
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Remark 1.4 The same remark holds as Remark 1.1 for the physical interpre- 
tation of f ll,15p . If we approach 

a(t) I inf Q(x, t) - p (t), 0(t) t sup Q(x, t) - p (t), 

in particular, then from (11.151) we obtain 

1 d 



-^\v\ 2 + u\u\ 2 j dx = 0, (1.11) 

which is equivalent to the classical energy equality (ll.ip for the Navier- Stokes 
equations. 

For the stationary solutions of the Navier-Stokes equations with v(x, t) = 
v(x),p(x,t) = p(x), and uj(x,t) = to(x), as an immediate corollary of The- 
orem 1.2 we obtain the following result, which is previously obtained by 
Galdi([3]). 

Corollary 1.1 Let N > 3. Suppose (v,p) is a stationary smooth solution of 
the Navier-Stokes equations satisfying: there exists a constant p such that 



lim {\v{x)\ + \p{x)-p \) = 0, (1.12) 



uniformly, and 



Then, v = 0. 
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Lf(M 3 )nL 6 (M 3 ) z/iV = 3, 
1)6 2JV (1-13) 
L^(R N ) ifN>A. 



Indeed, we first note that J^- > for > 4. Thus, if > 4, then 

3 jy 2N 

L—(R N )nL°°(R N ) L^(R N )nL°°(R N ) by the standard ^-interpolation, 
and the fact v E ^^(R^) is guaranteed by the hypothesis (I1.12p . Therefore, 
the stationary version of (11. lip together with the assumptions fll. 12j> and 
f l 1 . 1 3 1) imply 00 = 0. Combining this with div v — 0, we have v = Vh for 
some harmonic function h. Therefore the condition f)1.12p implies v = 0. 

Remark 1.5 For AT = 3, if we assume (I1.12p . and replace the assumption 
(I1.13P by J R3 |o;| 2 (ix < 00, it is still an open question whether v = or 
not. In this case we know by the maximum principle from (12. ip implies that 
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Q(x) < p for all x G M 3 . Therefore, by choosing a j. inf xgR 3 Q(x) — p , and 
(3 t in the stationary version of Theorem 1.1, we obtain 



\oo\ 2 dx = lim / \VQ\dS. (1.14) 

Therefore, we find that the desired Liouville type theorem holds if 



lim / \VQ\dS < / \u\ 2 dx. 

Ft° J{Q(x)= Po +/3} JR3 

In the case of smooth solutions to (NS,E) in a periodic domain (11.121) and 
(I1.13P are not necessary in view of the proof of Theorem 1.2. Therefore, as 
another corollary of Theorem 1.2 and its proof we obtain the following: 

Corollary 1.2 Let N > 2. Suppose (v,p) is a smooth solution of (NS,E) on 
T N x (0, T), where T N = /7* N is the periodic domain in R . We normalize 
the pressure by assuming J jN pdx = 0. Define the Bernoulli function Q as 
in U.ty) . Then, for all t G (0,T) and for almost all real valued functions 
a(t),f3(t) we have Then, for all real valued functions a(t) < (3(t) defined on 
(0,T) we have the following equalities. 



a(t)<Q(x,t)<p{t)} 



--M +„m )dx 



v \ \VQ\dS-v / \VQ\dS. 

>{Q(x,t))=a(t)} J {Q(x,t)=/3(t)} 



;i.i5) 



Remark 1.6 Similarly to Remark 1.4 the usual energy equality in follows 
immediately by approaching 

a(t) | inf Q(x,t), /3(t) t sup Q(x,t). 

We emphasize that we consider the energy equalities only for smooth solu- 
tions of the Navier-Stokes and the Euler equations. For the studies of the 
energy for the weak solutions of the Euler equations we refer to the results, 
among others, by De Lelilis-Szekelyhidi([2]), Shnirelman(j4]) and the refer- 
ences therein. 
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2 Proof of the Main Theorems 

Although the lemma below is well-known previously for the stationary solu- 
tions of the Navier-Stokes equations, we present its proof for reader's conve- 
nience. 

Lemma 2.1 If(v,p) solves (NS,E), then Q = p + \\v\ 2 satisfies 

~\v\ 2 + v\u>\ 2 = vAQ - v VQ. (2.1) 
Proof Multiplying (NS,E) by v, we obtain 
19, l2 1, I2 \ " N 



At; = (^l 2 ) -^(fy*) 2 - (2.2) 

^ ' .7,fc=l 



We compute 

AT jV AT 

j,k=l j,k=l j,k=l 
N 



2 



where we used the well-known formula Ap = —^f^idjVf-dkVj, which is 
obtained from (NS,E) by taking div(-). Plugging ( \2 3\\ into ( 12.2ft . we have 
(O). □ 



Proof of Theorem 1.1 We use argument similar to the one used in pp. 
We define Q(x,t) := Q(x,t) - Q Q (t). Then, 

lim Q(x,t) = Vt E (0,T) (2.4) 

uniformly. From (12. ip we have 

i/AQ-v- VQ = J^>| 2 + z/|u;| 2 . (2.5) 
2 or 

We assume below 

sup Q(x,t) > Pit) > a(t) > 0. 



The other case, > Pit) > a(t) > mi x€R N Q(x,t), is similar. Let us set 
D = {x E R N \a(t) < Q(x,t) < p(t)}. Thanks to the Sard theorem and 
the implicit function theorem <9D consists of smooth level surfaces in M. N 
except the values of a(t) and Pit), having the zero Lebesgue measure, which 
corresponds to the critical values of the function z = Q(x, t). It is understood 
that our values of ct(t), Pit) below avoid these exceptional ones. We integrate 
( 12. 5p on D, and use the divergence theorem to obtain 

' ^ 'v\ 2 + u\uj\ 2 1 dx = v I AQdx - I v ■ VQdx 



2 or 

= v I h-VQdS- I v-hQdS 
Jan Jan 

:= h + h, (2.6) 
where n is the outward unit normal vector on <9D. 

J {Q(x,t)=/3(t)} |VWi| J {Q(x,t)=a(t}) l V ^l 



V 



[ \VQ\dS-u [ \VQ\dS. (2.7) 

J(Q(x,t)=8(t\) JtQ(x,t)=a(t\) 



' {Q(x,t)=P(t}) J {Q(x,t)=a(t}) 

Next, by the condition divt> = we obtain 



n ■ vQdS — I h ■ vdS 

{Q(x,t)=P(t)} J {Q(x,t)=a(t)} 



= f3{t) \ (— n) ■ vdS — a{t) I n ■ vdS 

J {Q(x,t)=/3(t)} J{Q( x ,t)=a{t)} 

= (3(t) / divvdx — a(t) / divvdx = 0. 

J{Q(x,t)>P(t)} J{Q(x,t)>a(t)} 

(2.8) 

Combining <^M)-^M), we obtain (Oj) . □ 

Proof of Theorem 1.2 We set Q(x,t) := Q(x,t) — Po(t) here. Let 
a(t) < (3{t) be real valued functions on (0, T) with a(t) < Q(x,t) < Pit). 
The case a (£)/?(£) > is already covered by Theorem 1.1. Here we assume 
that 

m(t) := inf Q(x,t) < a{t) < < /3(t) < sup Q{x,t) := M(t). 
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Since 



{a(t)<Q(x,t)</3(t)} 

J 

J{a(t)<Q(x,t)<M(t)} 

f 

{P(t)<Q(x,t)<M(t)} 



1 9 , l2 i ... 

2 at 1 1 1 



- / ( \tt\v\ 2 + ^|w| 2 ) dx 

./{/9(t)<Q(x,t)<M(t)} \2ot J 

= [ (l^-\v\ 2 + u\uj\ 2 )dx+ I \VQ\dS 

J{a(t)<Q(x,t)<M(t)} K^Ct J J {Q(x,t)=/3(t)} 

by application of Theorem 1.1, it suffices to show that 

/ (lirM 2 + »W 2 )dx = - f \VQ\dS (2.9) 

J{a(t)<Q(x,t)} \^Ot J J {Q(x,t)=a(t)} 

for m(t) < a(t) < 0. The set D£ := {x G M N \ a(t) < Q(x,t)} is unbounded, 
while := {x E ~§t N \ a(t) > Q(x,t)} is a bounded set. Moreover, 

831 = dW U {oo}. 

Let R > be large enough so that dD^ C B R (0) := {x e R N \ \x\ < R}. We 
set := [Q + a] + = max{<2 + a, 0} below. We introduce the radial cut-off 
function a e C^(R N ) such that 



1 if \x\ < 1, 
if \x\ > 2, 



(7(|X|) = 

and < cr(x) < 1 for 1 < |x| < 2. Then, for each R > 0, we define 
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Multiplying (12 .5p by a R , and integrating on D?, we obtain by the divergence 
theorem, 



£ J^M 2 + v\u\* ) f*k<t-r ■ v I (T h >AQ+ ( ].r- / rr H r ■ Q+,].r 



v\ h-VQlo R dS -v VQl-V(j R dx- v-hQ+a R dx 
I Qtv-Va R dx 



:=h + I 2 + h + h, (2.10) 

where h = — is the outward unit normal vector to <9B°: . As in the proof 

I v y a | 

of Theorem 1.1, we have 

h = -u [ ^L.VQ* (7R dS = -v [ \VQ + a W R dS 
Jan 01 vyj 



-v 



[ \VQ\a R dS = -v [ \VQ\dS, 

J {Q(x,t)=a{t)} J {Q(x,t)=a(t)} 

(2.11) 

since a R = 1 on B R (0) D 831 = {x G R N \ Q(x,t) = a(t)}. By the fact 
<2q = on <9D" we have 

J 3 = 0. (2.12) 
Applying the divergence theorem again, we have 

J 2 = -v I VQ-Va R dx = -v Qh-Va R dS + v QAa R dx 

JO™ J dOt JD° 

— v \ QAa R dx, 



+ 



since Va R = on B R (0) D <9B". Hence, we have the estimates, 



\D 2 au L 



h\ < v \Q\\Aa R \dx < V — ||Q|| / dx 



R 2 11 n LN-2 (R<\x\<2R) 



{R<\x\<2R} 



N-2 

< Cv\\D 2 a\\ L o° ( / {\p-Po\ + -\v\ 2 )^dx 

{R<\x\<2R} 1 

< Cv[\\p-po\\n +\\v\\ 2 2n 1^0 

1 U LN~2( R <\ X \< 2 R) (R<\x\<2R) ' 
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as R — > oo, which follows from the assumption v G L°°(0,T; L 
and the Calderon-Zygmund inequality, \\p 



2N 
N-2 



where 



Poll. " < C|k|| 2 _aA!_, 

LN-2 L77-2 

we used the well-known velocity-pressure relation, p — p = Y^f,k^^' k ( v j Vk ) 
with the Riesz transforms Rj, j — 1, • • ■ , N in M. N . For the estimate of I4, 
since <9ID>° C Br(0), we have 



Q + a)v ■ Vo-rcIx 



{R<\x\<2R} 



< —\\Da\\ L o 



Qti ■ Vo~Rdx 



< c 

< C 




{R<\x\<2R} 



{R<\x\<2R} 



^HP-Pol + r|«| 



1 M L*r=T(.R<|a:|<2.R) 



3iV + \\V\\ 3N 

l J 2(N-i) (R<\x\<2R) 



\ V \\ -Mr \\ V \\ 3JV 

1 " (R<\x\<2R) " L^=T 



I ty 1 1 3iv 



L 7i r=T(.R<|:E|<2.R) 



-> 



JV 
JV-1 



dx 



(2.13) 



as R —7- 00 by the hypothesis lll.8p . where we used the fact a J{ R< \ x \ <2 r} v ' 
Va^dx = in the first line. On the other hand, if we pass R — > 00 in 
the left hand side of 1 12. 10p . we note that J oa \j^\v\ 2 a R dx — >■ J DQ |Jj|t>| 2 Gb 

by the dominated convergence theorem, while L a \u\ 2 aRdx — > L a |a;| 2 <ia; by 

the monotone convergence theorem. Therefore 112.91) follows from 1 12. 10p by 
passing R — > 00 in both sides of the equality. □ 
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